function [As bs] = support(A, b, w, h)
As = [];
bs = [];
[m n] = size(A);
for i = 1 : m
    [~,x,fval,exitflag] = evalc('linprog(zeros(1, 2), A, b, A(i, :), b(i), [0 0], [w h])');
    if exitflag == 1
        As = [As; A(i, :)];
        bs = [bs; b(i)];
    end
end
end
